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, solution of a random tiling model possessing an eight-fold symmet- 

ric phase. The solution is obtained using Bethe Ansatz and provides 
closed expressions for the entropy and phason elastic constants. Qual- 
itatively, this model has the same features as the square-triangle ran- 
dom tiling model. We use the method of P. Kalugin, who solved 
the Bethe Ansatz equations for the square-triangle tiling, which were 
found by M. Widom. 
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1 Introduction 



Random tiling models are ensembles of coverings of the plane, without gaps 
or overlaps, with a set of rigid building blocks or tiles. The most well- 
known application of such a model is perhaps the anti-ferromagnetic Ising 
model on the triangular lattice, of which the ground state configurations 
are described by rhombus coverings of the plane. 

The discovery of quasi-crystals® revived the modeling of structures by 
tiling models. Quasi-crystals are solid metaUic alloys without periodicity 
but with long range order (i.e. they have sharp diffraction peaks). Two 
competing scenarios have been put forward to explain the origins of such 
quasiperiodic ordering. One of them assumes a microscopic Hamiltonian of 
which the ground state configuration is quasi-periodic. This would be some- 
thing akin to the local matching rules used by Penrose. These deterministic 
tilings are perfectly quasiperiodic in the sense that their diffraction peaks 
have zero width. As an alternative there is the random tiling picture, in 
which one abandons the strict quasiperiodic long range order. These models 
do have matching rules in the sense that tiles may neither overlap nor have 
gaps between them, but these are not strong enough to enforce a unique 
behavior at large distances. This then gives rise to an ensemble of different 
tilings, or atom packings, inside a fixed volume. Assuming that all these con- 
figurations are (nearly) degenerate in energy such an ensemble has a nonzero 
entropy density. This entropy can reduce the free energy and lower it w.r.t 
the crystalline phase, which is supposed to be stable at zero temperature. 
This phenomenon thus provides a simple explanation of the stability of the 
quasi- crystalline phase. 

Of special interest in two dimensional random tiling models are those 
that fulfill the following constraints: (i) each kind of tile has a unique atomic 
decoration, (ii) each member of the random tiling ensemble is an allowed 
packing of atoms and (iii) each arrangement of atoms forming an equally 
good packing belongs to the tiling ensemble. These tiling models are called 
physical or atomistic. Examples of these are the binary tenfold rhombus^& 
and the twelve-fold square-triangle^iil' random tiling models. Recently Cock- 
ayne introduced a physical random tiling model with an eight-fold symmetric 
phase which he doped the octagonal analogue of the square-triangle tiling, 
despite of some distinctive characteristics. 13) In this paper we show that the 
comparison certainly holds concerning the exact solvability of the random 
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tiling model. We mean by this that, like for the square-triangle tiling,®EP 
it is possible to calculate the entropy exactly in part of the phase diagram. 
We obtain closed expressions for the residual entropy as well as for the pha- 
son elastic constants. The main results of this paper have already been 
published elsewhere, Q here we present the full calculation. 

The paper is organized as follows: In section ^ we present the tiling 
model and give a transfer matrix description to be able to enumerate all 
possible tilings. In section ^ we give a brief review of phason strain concepts. 
The results of this section will be used to identify our results and to obtain 
expressions for the phason elastic constants. Section ^ describes the diago- 
nalization of the transfer matrix by Bethe Ansatz. The resulting nonlinear 
coupled equations are solved in section |^. Finally we calculate the entropy 
density from this solution in section |^. 



2 The square-hexagon tihng and the transfer 
matrix 

The model that Cockayne described is a tiling of the plane by squares and 
various hexagons. By placing discs of radius 1/2 on each of the vertices, one 
obtains a physical packing. The hexagons can be thought of as being built out 
of two isosceles triangles and several rectangles depending on the length of 
the hexagon. The rectangles have sides 1 and a/2. Likewise, by drawing one 
of the two diagonals, the squares can be viewed as two triangles. The square- 
hexagon tiling thus is equivalent to the ensemble of vertex configurations of 
the rectangle-triangle tiling. The latter ensemble in turn is equivalent to the 
rectangle-triangle random tiling ensemble (see Fig. |1|) provided that the two 
ways of forming a square out of two triangles are counted as one.® 

Like for the square-triangle tiling we make use of the transfer matrix to 
calculate thermodynamic quantities like the entropy. The transfer matrix 
T is introduced by decomposing the tiling into layers. Different layers are 
bounded by the short horizontal edges, the horizontal diagonals of the squares 
and the almost horizontal diagonals of the ±7r/4 tilted rectangles. Here, ± 
refers to the tilt of the short edge of the rectangle w.r.t the horizontal axis. In 
addition, the layer boundaries cut the triangles and rectangles with a vertical 
long edge in half. The horizontal diagonals of the squares are denoted by the 
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dashed lines in Fig. |I]. In a similar fashion the tiling can be decomposed in 
'columns'. Different columns are separated by the short + edge, the vertical 
short edge and the almost vertical diagonal of the — tilted rectangle. The 
rectangles with a vertical long edge are cut by a column edge from the lower 
left to the upper right corner. In this way the tiles are deformed in such a 
way that the vertices of the tiling form a subset of those of the square lattice, 
see Fig. ||. 




Figure 1: Patch of the tiling model. The diagonals of the squares are not 
drawn. Dashed lines indicate parts of the layer boundaries. 

On this regular lattice the definition of the transfer matrix T is obvious: a 
matrix element Tij is if layer j can not be followed by layer i. Otherwise it 
is given by the statistical weight of layer i. In this paper we shall adopt the 
convention that T acts in the downward direction. 

In the following we will call the tilted rectangles R± and the rectangles 
with the short and long horizontal edge by Rg and Ri respectively. Because 
different tiles of the original tiling are mapped onto the same shapes on the 
square lattice, we have to decorate the new configurations. This is done with 
bold dashed and solid lines, see Fig. ||. Thus it is apparent that the horizontal 
short and long edges of the original tiling form domain walls between regions 
of i?+ , which we denote by type s and 1 respectively. If we think of the vertical 
coordinate as a time coordinate (increasing in the downward direction), the 
transfer matrix can be seen as an evolution operator for these domain walls. 
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The transfer matrix T is block diagonal with sectors parameterized by the 
numbers of domain walls, which are conserved under the action of T. 




Figure 2: Corresponding patch on the lattice. Bold solid lines are domain 
walls of horizontal short edges, referred to as s-walls in the text. Bold dashed 
lines represent the 1-walls of horizontal long edges. 

Between two layer edges on the square lattice, the s-walls step one unit to 
the right and the 1-walls do not move. When two walls cross, the s-wall 
may either jump over the 1-wall moving two places to the right and thereby 
creating a rectangle or, over two layers, the walls may exchange place 
creating a rectangle Rg. In the latter case the crossing therefore is completed 
after application of the transfer matrix twice. It may also happen that two 
walls of type s and one of type 1 cross simultaneously over two layers. The 
1-wall and the s-wall nearest to it then exchange place, while the second s- 
wall jumps over both these walls moving three places to the right creating a 
rectangle R-. 

We can express the tile densities in terms of the domain wall densities. 
We shall denote the horizontal size of the tiling by L and the corresponding 
system size of the lattice model by N. Let Ais = n^^ — n^j^ = 0, i.e. both 
types of collisions of two domain walls occur with the same frequency, and 
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let 2ns and ni be the number of s- and 1-walls. Apply the transfer matrix 
p = 2N — 2ns times on some initial configuration of domain walls at t = 
on the lattice, and suppose that both types of collisions occur for every pair 
of s- and 1-walls. The final state at t = 2N — 2ns will then be the same 
configuration of domain walls as the initial one shifted by —2ns. Counting 
the total number of tiles in this patch generated by T, we find the following 
relations between the wall densities and the rectangle densities per layer: 

p{nR^ + nn^ + 2nR_) = 2p{nR^ + nR_) = irisni. (1) 
A± = ur^ -nR_ = N - 2ns - ^i- (2) 

Similar expressions can be obtained for the different triangle densities. The 
total number of rectangles and of triangles per layer can then be calculated 
to be: 

rirect = N -2ns-ni + Ansm/p. , . 

ntri = 2(2ns + ni) - 12nsn\/p. 

The tile densities that belong to the quasi-crystalline phase are n^^ct/N = 
6 — 4a/2, ritri/N = 12^/2 — 16, corresponding to an area fraction of triangles 
at = 1/2. As a function of the domain wall densities, the model displays 
two incommensurate phases. A four-fold symmetric phase is formed in the 
high density region, at > 1/2, where the triangles form octagonal and square 
cells bounded by domain walls consisting of rectangles. There is a two- 
fold symmetric phase in the low density region where the rectangles from 
rectangular cells bounded by the domain walls consisting of triangles. 

As already mentioned above, the quantities ns and ni are conserved by the 
action of the transfer matrix T. To control the average value of Ais, the tiles 
Rs and Ri are given a weight exp(— 0) and exp(0) respectively. Furthermore, 
as the tiles Rs and i?_ in the lattice representation have an area that is twice 
that of the other two transformed rectangles, we have to introduce a chemical 
potential for them to compensate for this asymmetry. The tiles Rs and R- 
therefore get an extra weight exp(?7). 

The free energy per layer of the lattice model is given by the logarithm 
of the largest eigenvalue of T: 

/(ns,ni,0) = -log A 

= -NaN - 0Ais - r] [nR^ + nR_) , (4) 

where ajsi is the entropy per site of the lattice model. In section ^ will be 
shown that we can actually diagonalize T using coordinate Bethe Ansatz. 



6 



3 Phason strain 



In this section we will briefly review some of the phason strain concepts. 
Details can be found in a review on random tiling concepts by Henley. 1i2P. 
The position of every vertex in the tiling is of the form: 

3 

e|j' = (cos7rj'/4, sin7rj'/4). 

The tiling can therefore naturally be embedded in a four-dimensional hyper- 
cubic lattice using 

3 

(r,h) = ^n,(elU,^), .g. 
ej" = (cos 57rj/4, sin57rj/4). 

The vectors are chosen such that = (e|^',e^) are orthogonal. The con- 
tinuous height field, or representative surface, h(r) is as horizontal as possible 
for a perfect quasi-crystalline tiling. A generic member of the random tiling 
ensemble thus is described by fluctuations around this fiat surface. By inte- 
grating out the short-wavelength fluctuations one gets a smoothed function 
h(r), the phason field, that varies on length scales much larger than the tile 
edge length. The phason strain tensor then is defined by linearization: 

E = Vrh(r). (7) 

The macroscopic quasi- crystalline state corresponds to E = 0. The random 
tiling hypotheses now states that the entropy density cra(E) = S{Fi)/A has 
its maximum at E = and that it is quadratic in E near E = for any 
maximally random tiling. Like the square-triangle tiling, the octagonal tiling 
has the irrotational property. ^EJ) This implies that the most general quadratic 
form of the entropy density Ua consistent with 8-fold rotational symmetry is 

= ^a,o - (TrE)' + ]^K^ det E + O (E^) (8) 
We denote the deviations of the quasi- crystalline densities by 

A± = - ni - 2n,. 
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The quadratic forms in can be expressed in these, using L = n\\/2 + 2n^ 




The conditions on the elastic constants for E = to be a local maximum of 
a a are: 

ir^>0, irg>0, ^K,,-K^>{). (11) 

4 Bethe Ansatz 

It is possible to diagonalize the transfer matrix by making an Ansatz for its 
eigenvectors. Details of this so-called Bethe Ansatz are given in Appendix ^ 
A technical point we mention here is the fact that the s-walls split up in two 
different species, odd and even ones, both of which are conserved under 
the action of T. Their densities are given by Uo and respectively, with 
no + ne = 2ns. The eigenvectors and eigenvalue can be expressed in three 
sets of complex numbers corresponding to the three different kinds of domain 
wall. The eigenvalue is given by 

rio rtc 

A=n«^n^^' (12) 

i=\ j=i 

if the numbers Ui and vj satisfy the following equations 

k=l 

ni 

vf = {-r-'Wi^s+^'S'O- (13) 

k=l 

no no 

w;"" = (-)"'-^n(^'^«'+^'"Vv-')n(^s-+«'A"'^."')- 
i=i j=i 
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The equations (0) are the Bethe Ansatz equations (BAE). Using the sub- 
stitution 

^e,j = v], ^o,i = ul i;j = -w-'e^-'t (14) 
their logarithmic version can be written as 

N + rii 

^fsKoi) =0 (mod2jri), 

^^F.(U) = (mod2.i), (15) 

NFiiipj) = (mod27ri), 

where we have defined the following functions 

F,{z) = log(z)-^^J]log(z-^,)-^<^. (16) 

i=i 

X 1 '^'^ 

F,{z) = log(-^)--^log(eo,-^)-^$^log(ec,-^) 

1=1 i=l 

+ ^ log (A) + -<P - V- (17) 

The expression for the eigenvalue A of the transfer matrix (12) becomes 

no ric 

2 log (A) = log (^o,) + Yl log fe,) • (18) 

i=l i=l 

The Ansatz ( |60D together with the BAE (|1^) correctly give the eigenvector 
with eigenvalue (|18D of the transfer matrix for almost all densities. The 
description fails for the two states in which the lattice is filled completely 
with either s- or 1-walls. In these cases it is possible to have 'virtual' 1- or 
s-walls respectively, running around the cylinder. 



5 Thermodynamic limit 

It is clear that for the largest eigenvalue we must have Uo = rie- In that case 
the solution for ^o,i and ^e,i must be the same. If not, the largest eigenvalue 
would be degenerate which cannot be the case by the Perron-Frobenius the- 
orem. We therefore from now on omit the subscripts o and e and we set 
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TIq — Tig — Tig , cLS well as ^o,i = ^e,i = ^i- The roots and ijjj approach two 
curves in the complex plane (see Fig. ^. We shall refer to these curves as S 
and \E' respectively. 




Figure 3: BAE curves. 

From numerical calculations we note that there are no holes in the solution 
for the largest eigenvalue. Under the assumption that there are no holes in 
the thermodynamic limit of the solution we infer from the BAE ( [T5| ) that for 
{N oo) 

^^^/s(a)(a+i-a) = 2vri, (19) 
Nfiiijk){^k+i-i^k) = 27ri, (20) 

where fs and f\ are the derivatives of the functions Fg and Fi respectively. 
These equations allow us to convert the sums in (|TB|) and (|r^) into contour- 
integrals. By introducing the densities Qy = n\/N and Qs = Ug/N we obtain 
the following relations between /g and f\ in the thermodynamic limit: 

m - --^-T^d^, (21) 
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z 2711 J^, z-i 

The points &i, fej, 62 and 63 are the thermodynamic hmit points of V'l; V'np 
and ^1. From the form of equation (^) and ( p2D we see that /s(2;) is analytic 
along the contour Fg, though fi{z) is not. The function f\{z) differs on both 
sides of the curve S by the amount 

! + ^^''^--(l + QO/.W. (23) 



27ri Jjg_,,|=, 2:* - ^ 

where z^ is the point where Fg crosses the curve H. The analytic continuation 
of the function f\{z) along F^ therefore equals the function f\{z) + {l+Q\) fs{z) . 
A similar argument can be used for the analytic properties of the functions /g 
and /i along the contour Fi. The analytic properties of a linear combination 
G{z) = asfs{z) + a\f\{z) along the two contours are described by the so called 
monodromy operators Fg and Fi: 



a,\ 1 a. 



(25) 



If the two curves close, i.e. they have the same limit point 61 = 62 = b, it 
is not possible to close the contours Fg and F]. The only nontrivial closed 
curve is one that crosses both the curves H and \E', so that we are left with 
only one monodromy operator: 

r,r.= ( ^ ^7'). (26) 

V 1+Qi / 

which has the property that (FiFg)"^ = 1. This means that in this case the 
function G{z) is a single- valued function of the parameter 

t(,) = ( f'[-} )\ z{t) = b ^ (27) 
^ ^ \l-b* zj ^ ' l + bb*'^'^ ^ ' 

From now on we will write b = i|6|e''''. If we look at the orbit of fi{z) under 
FiFg we can see that fi{z), {1 + Qi) fs{z) - fi{z) , -fi{z) and -(1 + Qi)/g(z) + 
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fi{z) belong to different sheets of tlie Riemann surface of the same function 
G, where G is given by (Appendix jB]): 

G = ((e-'''/H + e'^/H-')-Qi{e'^/H + e-'^/H-'))^. (28) 



6 Entropy 

The knowledge of the function G allows us to calculate the eigenvalue from 
the BAE (|T^), see appendix |C|. Equations (^ and (^ show that the entropy 
per site for Ais = is given by: 

a^ = N'HogA-^^V- (29) 

This expression is precisely given by the solution ( ffSD of the BAE. 

At the QC-point, 7 = O,C = O,0 = O (for the definition of C, see 



appendix 0) and the domain wall densities given by ( p6|) and (67) become 

Qi = V2-l, 2Q, = 2-V2. (30) 
This implies via (|^) that the area per site A/N on the QC-point is given by 
A/N = I2V2 - 16, Qrect = 6 - 4^2, Qtri = 12^2 - 16. (31) 
Therefore the entropy per area at the QC-point is given by 



1 + V2 
0.1193642186 



a,,o = (A/iV)- V = i^^(log4-v^log(l + v^)) (32) 



Away from the QC-point, there are two incommensurate phases. There 
is a two-fold degenerate four-fold symmetric phase in which one type of tri- 
angle (say the up and down pointing ones) dominates while the other type 
is suppressed. The triangles form rectangular and octagonal cells bounded 
by domain walls consisting of rectangles. Each orientation of the rectangles 
occurs equally often. A four-fold degenerate two-fold symmetric phase is 
present, if one type of rectangle dominates (say as in Fig. |l|) forming 
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rectangular cells bounded by domain walls consisting of triangles. The 90 
degrees rotated rectangle (here is suppressed in this phase and the other 
two types of rectangle occur equally often, as well as both types of triangle. 

Two independent parameters describe the deviations from the QC-point 
into these phases. The four-fold phase is described by the parameter 7 which 
measures the difference in the densities of the two types of triangle. This 
phase also requires A^^ = 9 log A/90 = 0. The two-fold phase is described 
by e = A±/N = 1 — Qi — 2Qs. This parameter measures the difference in 
densities of the ± tilted rectangles, see equation (||). In this phase = and 
a log A/90 = 0. 

6.1 Two-fold symmetric phase 

In the two-fold phase 7 = Ais = = and the densities (|66|) and (|6^) 
become 

Q, = {V2-l){l+m{C) + 4'^{C)). 
2Q, = 2-V2-2V2^{C) + 2{2-V2)'i{C). (33) 

Using the fact that 5is/N = QiV2 — 2Qs = in the two-fold phase, C can be 
expressed in e: 

«(C) ^ ^(C) ^ (34) 

Since a nonzero value of e implies that C 7^ it thus not only affects the 
value of the residues, or equivalently the densities, but also the topology of 
the Riemann surface of G (see appendix P). The formula ( |64D gives the 
lowest order corrections compatible with the modifications to the residues 
using the above value for C. It does not however take into account the 
changes in the topology. Those corrections become comparable with (|5^) 
when 1 61 — 62!^''^ ~ e (see the end of appendix P). Since we want to know 
the corrections to the entropy density we need to know how small changes 
in the topology modify the integrals in appendix y. These corrections arise 
from the contributions of the poles at z = and z = 00 on other parts of the 
Riemann surface, which were previously disconnected. Since the net total 
residue on each part of the surface is zero, the lowest order correction to the 
integrals in appendix will be of the order of |6i — ~ e^- The entropy 
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density is linear in the integrals of appendix |C|, the lowest order corrections 
to the entropy can therefore be calculated using the modifications to the 
residues only. The entropy per area up to second order in e then becomes: 

= aa,o - e2i±}/? (^log4 + ^\og{l + v^)) . (35) 

Comparing this expansion with the expression for da in terms of the phason 
elasticity (|^) one obtains for the elastic constants: 

K^{y2-iy + K^V2= (V2-1) (log4 + V2 log (^1 + ^2) ) . (36) 
6.2 Four- fold symmetric phase 

In the four-fold symmetric phase e = 0. In this case the entropy throughout 
the entire four-fold phase can be calculated exactly. First we note that Ais = 
dA/d(f), which can be calculated analogously to the lowest order corrections 
in e. It is then found that 



OA dA 
dC 



= 0. (37) 

c=o 



The densities in this case are 



Qi = (^1 + ^2 cos I) ' (^1 - v^sin|j . (38) 

2Q, = (^l + y2cos|) ^ (^V2cos| + y2sin|) . (39) 
and the layer separation is 

A,N^l-iV-2-ir^p^. (40) 

v2 cos 7/2 + 1 

The area fraction of triangles in terms of 7 then can be written as 

<5tri| , /TT , ^,^2- cos 7/2 . 

at= . , = (v2 + 1) -i—, (41) 

A/N ^ ^ I + COS7/2 ' ^ ^ 
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so that the four-fold symmetric phase corresponds to 1/2 < at < 1. The 
entropy per area in this regime is found to have the following form 



2 + V2 



log 



4 



4cos^7/4 \ ' \cos7^ 
+ cos(^ + 2)logt<m(| + 2) 

/TT 7\ , /IT y 

+ cos(---jlogtan(--- 



(42) 
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Figure 4: Entropy as a function of a^. Sohd hue corresponds to the exact 
solution. Dots are numerical results for N=198. 



Expanding a^, up to lowest order in 7 yields 
2I + 



^a,0 - 7'^^ (4 - log4 - ^log (1 + V2)) . 



(43) 



To compare this with the phason strain description we note that in lowest 
order in 7 

7(^2 + 1) 



(44) 
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Using (|T0]) and (|^) we then find that 



4K^-K^ = ^^^^^ (^4-log4- V21og(l + V2) 

^ 0.4004597643... (45) 

Combining ([36| ) and ( ^Sj) one obtains the phason elastic constants: 



= 4(^2 - if ^ 0.2842712475 . . . (46) 

= ^^^?^(log4 + y21og(l + V2)-4(v^-l)4) (47) 
^ 0.7366252255... 

Qualitatively, this octagonal tiling displays the same behavior as the square- 
triangle tiling.® At at = I the model undergoes a so-called anomalous first 
order transition. This means that at fixed density the derivative of the en- 
tropy w.r.t. to the density, i.e. the pressure, has a jump, see Fig. ^. This then 
provides a mechanism for locking on the quasi-crystal phase over a range of 
pressures. 



7 Conclusion 

In this paper we successfully applied Bethe Ansatz to the square hexagon 
random tiling model. This model is a physical,® or atomistic, model of 
a quasi-crystal with eight-fold rotational symmetry. Amongst other things 
its solution answers the question posed in® if the random tiling hypothesis 
would hold for this model. It remains unclear however, what the precise 
conditions for solvability of these tiling models are. In this respect, the 
irrotationality constraint, which is satisfied by both the square-triangle and 
the square- hexagon tiling, might be an important guide. Such a constraint, 
not necessarily irrotationality, which reduces the degrees of freedom probably 
has some infiuence on the integrability. The tenfold tiling by rectangles 
and triang les© however does admit a Bethe Ansatz. We know not of any 
constraint for this tiling, nor have we solved the Bethe Ansatz equations. In 
the future we hope to present some numerical results for this tiling. 
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A Bet he Ansatz 

In this section we show how to set up an Ansatz for the eigenvectors of the 
transfer matrix T to solve the eigenvalue problem. The conserved quantities 
make T block-diagonal. We can make use of this fact by diagonahzing T in 
each sector separately. It is helpful to start with the easiest sectors 2ns — 0, 1, 
ni — 0,1. Later we generalize the method to arbitrary values of 2ns, ni. 

A.l Domain walls of one kind only 

The sectors where either ni = or ng = arc very simple. E.g., in a sector 
with ni = 0, T just shifts all s-walls one step to the right. If we assign a 
wave-vector pi to every s-wall and denote the coordinate of the ith s-wall by 
Xi, the eigenvectors ip{{xi}) of T in these sectors can be written as a product 
of plane waves: 

2ns 

^{{x^})=lle-'P^^^. (48) 

i=l 

The eigenvalue in terms of the wave- vectors is: 

A = ]^e'f\ (49) 

i=l 

When using periodic boundary conditions the wave-vectors have to obey the 
following equation: 

e'^P' = 1. (50) 

These equations are easily solved to give the spectrum of T in these sectors. 
For later convenience we will fix some notation. The wave-vectors for the 
1-walls will be called qj and in the following we will make frequently use of 
the exponentiated momenta Ui = exp(ipi) and Wj = exp(igj). 
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A. 2 The sector n\ = 1, 2ns = !• 



The simplest non-trivial sector is that with one s- and one 1-wall, see Fig. ^ 
Let z be the coordinate of the 1-wall, the eigenvalue equations then read: 



Aipi^x; z) = iIj{x — 1; z), x > z + 1. 

Kip{x; z) = ip^x — 1; z), X < z. 

AV^(x + l;x) = e^'^-'^'^/^'ilj{x;x)+e'^^{x-l;x). 

A^(x;x) = e^'^-^^/^'ilj{x;x + l). 



(51) 
(52) 
(53) 
(54) 



We do not have to specify exactly the positions of the walls in the interme- 
diate state ipix; x) when the walls cross via Rg. We only denote this state by 
■ip{x;x) which we can express in ip{x; x + 1) with equation (l5^). 





X X -\- 1 



x x -\- 1 



Figure 5: Crossings of two domain walls. 



We can solve this eigenvalue problem by making the following plane wave 
Ansatz for the eigenfunction 

^p{x;z) = A{p,'n)u~'''w-\ x>z. 

ip{x;z) = A{-K, p)u:;^''w~'' , X < z. (55) 

From this Ansatz, in which we have introduced the auxiliary labels vr and p 
for later convenience, follow the eigenvalue and the relation between the two 
amplitudes before and after the collision 

A = M^, A(p,7r) = S'is(p,7r)A(7r,p). 
5is(p, vr) = {e^u^ + e'^-^u-^w;^) . (56) 
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If we take periodic boundary conditions we get the following equations for 
the complex numbers u^, and Wp-. 



A{p,7l 



A{n,p) 

{u^Wpf = 1. (57) 
One can check that for arbitrary rii and 2ns = 1, the eigenvector is given 

by 

ni 

^(x;;2i,...,;2„J = J]A(r)M;^J]«;;f , A = u^. (58) 

p k=l 

Here we have introduced a shorthand notation F for the dependency of the 
amplitude on the permutation p and the order of types of domain walls. This 
information is coded in the following way. Let r be the vector of coordinates 
X, Zk of all domain walls, ordered so that < r^+i. The entries of the vector 
r are the elements of the permutation p together with the auxiliary label vr. 
The order of succession in F of vr and elements taken from p matches that of 
X and the z^s in r, like in equation ( |55D 

Amplitudes differing in an exchange of neighboring s- and 1-walls or an 
exchange of two neighboring 1-walls are related by the scattering matrix el- 



ements 5'is(pfe,7r) (|56D and Sii{pk, pi) respectively, where Sii{pk, pi) = —1 for 
k ^ I. Imposing periodic boundary conditions gives us the BAE for these 
sectors: 



u. 



N ^(Pl,---,Pni,7r) 



A{7l,pi, . . . ,p„J 



k=l 

N 



Sisik,n) H Snik,j) = i-r-'SUk,7r). (59) 
j^k=i 



A. 3 Sectors with arbitrary values of rig and ni 

Let z be the coordinate of the 1-wall and Xi the coordinates of the ith s-walls. 
By letting t be the vertical coordinate and q the number of 1-walls to the left 
of the ith s-wall, we see that {xi + t + q mod 2) is a conserved quantity for 
every s-wall. This means that the s-walls split up into two kinds: odd and 
even ones. We shall therefore change our notation a little bit: Xi and yj will 
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denote the coordinate of the ith odd and the jth even s-wall respectively. 
Their total numbers will be denoted by rio and rig and 2ns = ""-0 + ^6 and we 
shall use the shorthand 5*0,6 for S's„,so- 

The Ansatz for the eigenvector in a general sector is given by: 

Ipi^Xi, . . . , Xn^, yi, . . . , Uriaj ^1; • • • 5 ^ni) — 
rio no ni 

E E E ^(n n n n «^p7n m 

IT fj. p i=l j=l k=l 

with eigenvalue 

A=n^^n^^- (61) 

i=i j=i 

It turns out that in higher sectors all relations factorize into two-domain wall 
scattering relations, like (|56|). The BAE for Ui, Vj and Wk using periodic 
boundary conditions then are 

= n 'S'o,o(^,m) J^So.ol^, j) n'S'io(A;,i). 

m^i=l j=l k=l 

rio rio rii 

vf = HSoA^^j) n Se,ei3,m)l[SUk,j). (62) 

i=l m^j=l k=l 

rio "-c 'il 

Wk"" = \{Sio{Ki)\{S,,{k,j) J] Su(A;,m). 
Here, 5*10 is given by (^) as well as 5*10 with u replaced by v. Furthermore, 

•^ll = S'0,0 = •S'e.e = —1, 5'o,e = 1- (63) 

B Calculation of the function G 

In the following we will always assume that the two curves have the same 
limit point. In this case the Riemann surface of the function G becomes 
disconnected and consists of an infinite number of discrete parts comprising 
four sheets each. Each of these parts can be mapped onto the plane by (|27|). 
We fix a particular part and sheet of the Riemann surface by choosing G 
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equal to /i at 2; = 0, on the sheet with the point t = e™^^. The form of G on 
the other sheets of this part of the Riemann surface is then determined by 
the monodromy operators Fg and F]. The finiteness of the densities Qs and 
Q\ imphes that the forms fs{z)dz and fi{z)dz are non-divergent at z = b and 
z = b*. The only singularities therefore are the simple poles at z = and 
z = 00. The residues of Gdz at these points can be easily calculated from /i 
and its monodromy using equation (^Gf) . They are listed in Table |l|. From 
this table we see that the functions /i, (l + QO/s — /i, — /i and — (l + (5i)/s + /i 
belong to the different sheets of the same part of the Riemann surface of the 
function G. 



z 




G 









/l 


1 


00 


ie-i7/2 


(i + gi)/s + /i 







_g~i./4 


(i + gi)/,-/i 


Qi 


00 


_e-i7/2 


-/i 


i-2gs 





_gi^/4 


-/i 


-1 


00 


_ie-i7/2 


-(i + gi)/,-/i 









-(l + Qi)/s + /i 


-Qi 


00 


e-i7/2 


/i 


2Qs-l 



Table 1: Poles and residues of Gdz 
The form Gdz is defined unambiguously by its residues and equals 

8 

Gdz = y^^dt 

n=l ^ 

Y 2 2 

+C {t + + c* + t3) J-^dt. 
^ ^ ^ ^ / dt 

where C given by 



(64) 
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+ (l-2Q,)(l-i)e^^/2 + 2isin 



(65) 



Or equivalently, one can write the domain wall densities as functions of 7 
and C: 



2Qs 



Qi 



- 2^{C) (^cos I + cos + 253(C) (^sin ^ 
(^l + V2cos|) ^(^1-V2sin| 



. 7 
sm — 

2 



(66) 



+ 23fJ(C) (cos ^ + cos ^ + sin ^ 
y 2 2 2 



sm ■ 



, ^c-^r-^ / ^7 7 .87 

+ 2^(G ) ( cos — + cos — — sm — + sm 

Zi Zi 



(67) 



The monodromy properties of the curves S and \& only hold when fei = 62- 
These curves are defined by the form Gdz taking only purely imaginary 
values on the tangential vectors of H and see equations (^) and (pUj). 
Given the initial conditions these solutions must be unique except precisely 
at the points t = and t = 00. At these points the different solutions should 
meet. If this were not the case the curves would intersect or would not close, 
destroying the assumed monodromy assumption. The condition 3?(Gdz) = 
must therefore be satisfied identically at the points t = and t = 00. From 
(HD and (13) it follows that 



G^dt 
dt 

Opt 

dt 



4C(1 + e^''')dt 



t=o 



t=oo 



t=0 



-4C*(l + e 



-2i7 



dt 



(68) 



i=oo 



The monodromy assumption translates itself therefore in the condition C = 
0. The Riemann surface of G then has infinitely many disconnected parts 
with four sheets each. If bi 7^ 62 these parts become joined by regions of size 
\bi — b2\. One can still use the variable t with b = {bi + &2)/2. The function 
G then is single- valued except for the regions around t = and t = 00 where 
\t\ ~ 1 61 — 62 T'^^ or \t\~^ ~ l&i — ^2!^''^ respectively. Equation (|6^) gives the 
corrections to the changes in the values of the residues if C 7^ 0. It does 
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not however take into account contributions arising from the changes in the 
topology of the Riemann surface. For small t the form Gdz is approximated 
by terms of the order of t^dt and Cdt. As a lowest order approximation (|5^) 
thus is valid if ~ \bi - h2\^'^ < C. 

C Calculation of A 

In order to remove singularities at z = and 2; = oo occurring in some of 
the expressions in this appendix, we introduce the following forms: 

dt 

dt 

(69) 

dt 







r4 


+ 








-u 


t 


-ts 


gsdz = 










re 




-t2 


+ 


t 


-u 


hidz = 




n 








(r 


-ti 


+ 


t 


~h 


hsdz = 




rs 


+ 






(r 




t 





The BAE ( |T5| ) allow us to evaluate the following expression using the defini- 
tion Fi as given by equation ([T7|). 



POD 

Ji,oc = ^ (/i - 91) dz 
Jb 

= 3? lim {Fiz{t)) - r4 log(t - ^4) - rg log(t - tg)) 

= (1 - 2a) log ' + log A + 2(1 - Q,)(P - f]. (70) 

The equality sign on the second line follows from the fact that Fi{z) is a 
primitive function of fi{z) and that 6 is a solution of the BAE (|T5|). The 
expression on the third line follows from the definition of ^^1(2:) as given by 
equation (|l^). Similarly we find using also the definition for Fs{z), equation 

POO 

Js,oo = ^ / iil + Qi)fs + fi-9s)dz 
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(2 - Q, - 2Qs) log + AT-i log A 



+ 2(1 - Qi - gs)0 - r]. (71) 



1,0 



3fJ / (/i-/ii)dz 
4|6| 



log 



COS' 




N-HogA + 2{l-Q,)(f)-r]. 



(72) 



3f^ / ((l + Qi)/s-/i-/is)dz 
4\b\ 



Q\ log 



cos 7 1 



+ Ar-ilogA-2si-2(l + Qi-Qs)</' + ??. (73) 



with 



ni 



(74) 



si = AT-i^log^,. 

i=i 

Eliminating and si from the equations (l70|) — ([7^) we find for the eigen 
value A and the chemical potential t] 



N Mog A - 



QsQi 



QsQ\ 



■In 



QsQi 



-^1,0 - 1 



l,oo 



1-Qs 



log^— . (75) 



On the other hand, the integrals Ji^oo and Js^oo can easily be integrated in 
the t-plane using the results of appendix |B[ 



Jl,oc = '^nlogltn - tsl • 

8 

-^s,oo = ^ rn\og\tn - ^2 | • 



(76) 



(77) 



With the help of Table |T] of appendix |B| these expressions are given explicitly 

by 



l,oo 



- log 

2 ^ 



1 - cos (7r/4 + 7/2) 



1 + cos(7r/4 + 7/2) 



— log 

2 ^ 



1 + cos (7r/4 - 7/2) 



1 - cos (7r/4 - 7/2) 



.(78) 
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J. 



1 - COS (7r/4 - 7/2) 



l + cos(7r/4-7/2) 
And similarly J\ q and Jg^ give: 



1 - cos (7r/4 + 7/2) 



1 + cos(7r/4 + 7/2) 



.(79) 



1,0 



(2 - Qi - 2Qs) log 



1 + cos (7r/4-7/2) 



1 - cos (7r/4 - 7/2) 

1 + cos(7r/4 + 7/2) 



+ -(1- 2a) log 



J. 



s,0 



-(2-Qi-2g,)log 



1 - cos (7r/4 + 7/2) 
1 + cos (7r/4 + 7/2) 



1 - cos {tt/A + 7/2) 

1 - cos (7r/4 - 7/2) 



+ -(l-2Qs)log 



1 + cos(7r/4 - 7/2) 



(80) 
(81) 

(82) 



Combining (fTSI) — (P^D and (|7Sp then gives the eigenvalue expressed in the 
densities (given by (^) and (|67D), the chemical potential 77 and the parameter 

7- 
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